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ABSTRACT 

We investigate the feasibility of carrying out likelihood analysis on the velocities of 
galaxy clusters to determine power spectrum parameters. Using N-body simulations 
of cosmological density fields we show that the velocity field traced by clusters is 
highly non-linear on small scales (r < 10 h~^Mpc), mostly due to strong infall 
into superclusters. The transverse and line of sight velocity correlation functions 
relevant to likelihood analysis deviate significantly from the linear theory values. 
Unfortunately, the sparseness of cluster samples means that they are not amenable 
to any rigorous smoothing procedure to rectify this. We are forced instead to 
remove close pairs of clusters which would contribute disproportionately to the 
likelihood. This reduces the size still further of the small samples of observed 
clusters available to us, forcing us into the conclusion that using straightfoward 
linear theory covariance analysis on clusters is unpromising. Nevertheless we apply 
the technique to two sets of observational data and find marginal evidence for more 
power in the mass density field than standard CDM. The normalisation is consistent 
with Vf'^CTg, ~ 0.7 for both data sets, within the (large) errors. 

1. Introduction 

In linear theory, two-point statistics can be used to define the properties ol the 
velocity field entirely. If it were possible to find tracers of cosmological velocity fields 
that obey this criterion, then covariance analysis could be used to remove uncertainties 
due to sample geometry^'^. Clusters are potentially good tracers because velocities of 
galaxies in clusters can be averaged to give smaller errors. 

2. Comparing simulations of cluster formation with linear theory. 

2.1. Catalogues of clusters from N-Body models. 

A series of simulations were carried out in order to assess how well clusters trace 
the linear velocity field. A P"^M N-body code"^ was used to follow the evolution of 
cold dark matter and clusters were selected from the simulations by a percolation 
technique. This is as descibed in [4], except that here we have simulated a box size 
of size 600 h~^Mpc (where IIo=100h kms~^Mpc~^) with 4 X 10^ particles. The initial 
power spectra (hereafter PS) are those of [5], with the shape parametrised by F, 
equal to Oq^ in CDM models. We will present results based on 2 simulations of the 
Standard CDM scenario ( for SCDM, Oq = l,r = 0.5) and 4 simulations of Low 
Density CDM (LCDM has Oq = 0.2, F = 0.2, and a cosmological constant A such 



that A/3Hq = 1 — Oo)- The models are normahsed so that the rms mass fluctuations 
in 8 h~^ Mpc spheres, erg = 1.0 at the present day. Because we choose to study 
velocities of rich systems, a lower mass bound is applied so that the clusters have a 
mean separation of 30 h~^Mpc. 

2.2. Velocity correlations. 

Although the rms velocities of simulated clusters agree roughly with linear theory"*, 
approximately 5% of clusters occupy a signiflcant non-Maxwellian tail. Examination 
of the position of these objects in space shows many of them to be infalling into 
superclusters. To flnd the range of scales over which these non-linear effects may 
be important we examine the velocity correlation tensor, ^^^(r). This describes the 
covariance of different velocity components as a function of distance between two 
points. It can be written as a sum of the radial (^||(r)) and transverse ±{r)) 
velocity correlation functions^ : 

%(r) = (K-(x)y,(x + r)) = *i(r)% + [*||(r) - *i(r)]f,r,. (1) 

In linear pertubation theory these velocity correlations are given by integrals over 
the PS^: 

*||(r) = I -j^eM-k'R'j){jo{kr) - ^^)dK (2) 

*i(r) = ^^M-k'R}r-^)dk, (3) 

where j„ are the usual spherical Bessel functions, and the velocity held has been 
flltered with a Gaussian window of radius Rj. We will use Rf = 3h~*mpc in our 
comparison of linear theory with the velocity correlations of simulated clusters, as 
this value provided the best flt to the one-dimensional velocity distribution'*. 

We have evaluated ^||(r) and ^ i_{r) for samples of simulated clusters, and the 
results are shown in Fig. 1, together with the linear theory predictions. It can be seen 
that cluster pairs with separations < l{)h~^ Mpc show a strong tendency to fall toward 
each other (^||(r) is stongly negative) rather than streaming together as expected in 
linear theory. 

2.3. Maximum likelihood analysis. 

The relative probablility of observing a set of N line-of-sight peculiar velocities Vi 
at positions r; in a Gaussian held is given by a multivariate Gaussian distribution*: 

fl^V 1 
nyi)d''V = , ^ ,,,, exp(--K-^-V;0, (4) 

where Ci^ = {ViV j) is the covariance matrix of the sample. In linear theory, the terms 
in are simple functions of the velocity correlations deflned in equations 2 and 3. 
If clusters 1 and 2 are positioned at ri and T2 relative to the observer, and the angles 
61 and 62 are deflned by cos ^1 = ri • ri2 and cos 62 = ^2 • f 12 then 




Figure 1: The transverse (upper lines) and parallel (lower lines) velocity correlation 
functions for rich clusters in (a) the Standard CDM model and (b) Low Density CDM. 
The smooth curves are the linear theory prediction (equations 2 and 3) and the points 
are evaluated from samples of simulated clusters. 



Ci2 = cos 6'icos6'2*||(ri2) + sin 6'iszn6'2*i(ri2). (5) 

For a real sample with velocity errors, there will also be additional contributions 
to the covariance matrix. We will assume that these errors are uncorrelated and have 
a normal distribution, so that only the Cu terms are affected, each being increased 
by crf^ the estimated error on the velocity of cluster i. 

The PS used in the calculation of ^|| and ^ ±{r) is again that of [5]. We maximise 
the likelihood given by equation 4 in the parameter space of F and effective amplitude 
0°'^(T8. Whilst it would be desirable to test a fit to the PS which is a sum of step 
functions, this way of finding a more general answer is not practicable with this 
method, which involves many inversions of the matrix C . We are also limited by the 
small number of data points. The parameter F should be a fair phenomenological fit 
to the PS of most reasonable models, though, and we would expect its value to lie 
somewhere around 0.2, that inferred from the shape of the galaxy PS'''^. 

Testing the procedure on samples of linear theory velocity fields reveals that al- 
though there is rather more scatter in the shape parameter than amplitude, the 
method recovers the right values on average. We then carry out likelihood analysis 
on samples of simulated clusters. Disappointingly, we find the method unable recover 
the correct parameters. Tests have been carried out^ which show that clusters in 
close pairs dominate the likelihood and skew the parameter determinations to wildly 
incorrect values. Some means of mitigating the non-linear effects seen on small scales 
needs to be applied. Unfortunately we find that clumping the data can introduce 
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Figure 2: The distribution of best fit values in T — erg space found by likelihood 
analysis of 200 samples of simulated cluster velocity fields (60 clusters each in fOO 
h-^Mpc spheres). Panel (a) is for SCDM (F = 0.5,O°-V8 = 1) and (b) is LCDM 
(r = 0.2,O°'^(T8 = 0.38). Clusters with a nearest neighbour within fO h~^Mpc have 
been removed from the sample, as well as the fastest 3 per cent of clusters. 

large biases^. By considering grouped data as single objects we are removing small 
scale power from the velocity field. Unless we can smooth rigorously with a filter 
that can be introduced into the theoretical covariances of equation 5, we will recover 
values of F and erg which are systematically and unpredictably low. We therefore 
choose to excise clusters with close neighbours (r < 10h~^ Mpc) from the calculation 
altogether. As a final crude attempt to remove non-linearities, we remove the fastest 
3% of clusters on the basis of their position in the one-dimensional velocity distribu- 
tion. The distribution of best fit PS parameters for tests under these conditions is 
shown in Fig. 2. It is evident that some biases are still present which should be borne 
in mind later. We have also tried adding simulated observational errors. The results 
(which are not shown) show that scatter is increased but do not bias the results, as 
long as we use the real space separations of clusters. 

3. Applying maxmimum likelihood analysis to observed cluster velocities. 

3.1. The cluster samples. 

We deal here with two separate samples of observed clusters, one of which (AFFHM) 
is a compilation of 65 clusters^°'^^'^^'^"^'^'* and is described in [4]. The other (FP) con- 
sists of distance and velocity measurements of all 119 Abell clusters within 150 h~^mpc 
using the Brightest Cluster Galaxy method^^. These distances have fairly large (16%) 
errors, as only one measurement is possible per cluster. On the other hand, the large 
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Figure 3: Contour plots resulting from likelihood analysis of the cluster velocities in 
(a) the combined AFHMM sample and (b) the Lauer & Postman sample. The heavy 
contours enclose 68 per cent, 95 per cent and 99.7 per cent of the probablity if the 
distribution of — 2 In(likelihood) follows a distribution. 

volume occupied by the sample has already been used to infer a bulk flow which has 
significant power to rule out models^^'^'''^. 

3.2. Results. 

Applying the same treatment used on the simulated clusters in Sec. 2 to these data 
sets, we are left with only 38 and 89 objects in the AFFHM and LP sets respectively. 
We then calculate the likelihoods resulting from equation 0.6 on a dense grid of F 
and 0°'^(T8 values. The results are shown in Fig. 3 with the heavy lines representing 
confidence limits that the result lies within those contours. Both samples point to 
a best fit amplitude of O^'^CTs ~ 0.7 and to a PS with more large scale power than 
SCDM (F ~ O.f for AFHHM and ~ 0.05 for LP). If we assume that the PS shape 
is known to be 0.2 as derived from measurements of the galaxy distribution, then 
maximising the likelihood for variations in only 0°' ^8 gives 0°-V8 = 0.73tEJi? for 
AFHHM and n^-'^ag = 1.44^°;^^ for L&P (90% confidence hmits). 

4. Conclusions 

We find that simulated clusters show strong non-linear infall into high density re- 
gions. Clusters are distributed too sparsely to allow for smoothing with a well defined 
filter in an attempt to reduce the non-linear effects. Any crude grouping of the data 
will introduce strong biases into the determination of power spectrum parameters. 
Rather than abandon linear theory covariance analysis altogether, we have chosen to 
remove high peculiar velocity clusters from the likelihood calculation. Applying the 



technique to two different sampies of observed ciusters resuits in determinations of F 
and (78 wliicli are consistent (F ~ 0.05 — 0.3 and erg ~ 0.4 — f .0). Taking into con- 
sideration possible remaining systematic errors, botli tlie SCDM and FCDM models 
appear to lie ~ f sigma away from the best fit power spectrum, which itself is near 
to that of the Mixed Dark Matter model^®. 
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